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One dimensional Dirac equation is analysed with regard to the existence of exact (or closed-form) 
solutions for polynomial potentials. The notion of Liouvillian functions is used to define solvability, 
and it is shown that except for the linear potentials the equation in question is not solvable. 



- - . I. INTRODUCTION 

^ , 

■ The question of solving differential equations in explicit terms has lost some importance with the development 
Cs| ] of numerical methods, and it is also clear that solvable cases are non-generic in the class of all physical models. 

(— ( . Although the linear ordinary differential equations are mostly dealt with by series expansion or even used to define 
' new functions, it is still instructive if an explicit solution can be found. Obviously, exact solutions provide deeper 
^"5 , quantitative insight into whole classes of solutions (depending on some parameters); and in case of perturbative 

■ approach, where the linear equation is just the first approximation, they allow to shift the numerical approach to the 
' next order. Examples include Schrodinger like equations [H, f2|, and perturbation equations in cosmology fsj. 

When solving the Dirac equation it usually boils down to solving a second order linear differential equation - 
be it by an appropriate change of variables, or by a particular ansatz. Further assumptions or symmetries lead to 
I i special solutions which can be obtained explicitly, most of which are presented in Depending on the physical 
1^ ' theory one starts with, even though the form of the equations stays the same, the results on solvability differ. This 
happens because of various definitions of the potential - for example, in some supersymmetric theories the potential 
^ W — V"^ + V is taken to be a generic polynomial, whereas here we take V to be generic, which leads to restrictions 
on U. For results on integrability of such equations see for instance 

There are also many definitions of "explicit solutions" depending on which special functions are considered to be 
] simple enough. Here we will be using the Liouvillian extensions to construct solutions. Since the equation will be of 
. the form 

f"{x)^r{x)f{x), (1) 

CO ' where r{x) is rational, the sought solutions will lie in some extension of the field of rational functions over C. The 
, extension will be called Liouvillian if it is composed of a finite number of the following steps: 

, 1. Adjoining an element whose derivative is in the field. 
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2. Adjoining an element whose logarithmic derivative is in the field. 

3. Adjoining an element algebraic in the field. 

These functions include all the elementary ones (polynomials, exponential, logarithm), and also special cases of 
transcendental functions. The reason for introducing such class is that it naturally follows from the differential Galois 
theory Q and that for the above equation there is an algorithmic approach to checking for solvability @. To put it 
generally, when the solutions are Liouvillian, the differential Galois group is solvable, and there exist invariants of the 
equation for which one can look. This means that, at least in theory, it is possible to check a given class of equations 
and rule out the existence of any additional exact solutions - an advantage usually not found in other approaches. 

Thus, to put it simply, the aim of this letter is to find Liouvillian solutions of the Dirac equation into which a 
polynomial potential has been incorporated. As we will see, the conditions for such solvability are almost never met. 

II. DIRAC EQUATION WITH A POTENTIAL 

We take the one-dimensional form of the Dirac equation 

idtip — {ap + (3ni)^p, (2) 
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where p is the momentum conjugate to the coordinate x, i.e. p = —idx, and the Dirac matrices satisfy = 0^ — 1, 
{a, /?} — 0. This form follows from the standard four-dimensional form, when the bispinor ip depends on t and x only. 
It can then be taken as a two component vector ■0 = (^i (i, x), ■02(i, 2;)) for simplicity. Also the time derivative can 
be eliminated by the ansatz 'ip{t,x) = exp{—iEt)ijj{x) , so that the equation is 

Etp = {-iadx + I5m)ij. (3) 

The inclusion of the potential which is specified by only one scalar function can be performed in two ways. The 
polynomial potential V — Xx"' + . . . + Ao can be the time component of a four- vector - like the Coulomb part of the 
electromagnetic field. The equation would then become 

(y + E)i^ = {-iadx + Pm)^. (4) 

Alternatively, one could introduce scalar coupling, which modifies the m term to 

E^P ^ {-iadx + P{V + m))^. (5) 



In the first case, the explicit equations are 



V'a = -im + V)^J2 + Eijji, 



or, as second order equations. 



ijj'l = ([/' + C/2 - E^)ib 



i;'^ = i-U' + ~ E^)tP2, 
where U = m + V and the Dirac matrices were taken to be 



^-i / ' ^ U 
The vector coupling, on the other hand, gives 

provided that 



(6) 



(7) 



(9) 



and U =: V + E. It is obvious that the vector coupling can be transformed into the scalar one with 

V — >■ —iV, E —im, m lE. (11) 



III. THE SOLUTIONS 



Let us start with the scalar coupling. Since the spinor components are connected through the first order equations 
^ it suffices to check just the first of equations ([7]). For polynomial V , and thus polynomial U it is straightforward to 
apply the aforementioned Kovacic algorithm. Because the only coefficient in the equation is an even degree polynomial, 
only the first case to be considered in depth, and the possible solution will be of the form Pexp J wdx, with a monic 
polynomial P and a rational function lo. 

Since the only singular points is the infinity one needs to obtain the expansion [Vrjoo of r — U' + — E'^ . Because 
only the polynomial part of the expansion matters, it is fully and uniquely determined by comparing the coefiicients 
of r and [y/r^"^. Obviously, there are at most n -f 1 terms of the expansion, and the term alone fixes the n + 1 
highest terms of r. Thus, U itself is the solution and 



r - [VtooY = r-U^^ nXx"-^ + . . . - E^ , 



(12) 
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If n 7^ 1, the algorithm gives uj = U, and deg(P) = with the additional condition that E^P = 0. In other words, for 
E = Q the solution is 



If n = 1, there is only the zeroth order term in (fT2)) and for i?^/(2A) e N we have a whole family of solutions. This 
is the known case of the so-called Dirac oscillator solvable by Hermite polynomials [1]. 

For non-zero values of E no other (non-constant) polynomial potential gives Liouvillian solutions. Thanks to the 
relations (jlip . we can see that the same holds for the vector coupling, with the appropriate interchange of E and m. 
Namely, for a linear potential one has the Hermite solutions as above, and for m = the solution is 



The above results can be summarised as follows: 

Theorem 1 The one- dimensional Dirac equation with a polynomial potential of degree n > I only has Liouvillian 
solutions when E = (scalar coupling ^) or m = (vector coupling given by formulae (|13p and (|14p respec- 

tively. 

For n = 1 the solutions are expressible by Hermite polynomials given in f^], and n = is equivalent to the standard 
Dirac equation. 

Two remarks are in order. First, note that in the special cases m = {E — 0), the solutions given above hold for all 
potentials, not just polynomial ones as can be checked by direct differentiation. Second, the practical consequences 
of the outcome are that the spinor is not expressible as a polynomial function of position, as is often the case in 
solvable models of quantum mechanics. This still leaves the possibility of using new transcendental functions or series 
expansion approach, but precludes the direct (explicit) computation of norms or other general expressions involving 



This paper was supported by grant No. N N202 2126 33 of Ministry of Science and Higher Education of Poland. 



[1] Acosta-Humanez P. B., PhD Thesis, Universitat Politecnica de Catalunya f2009'). larXiv:0906.353"2l 

[2] Acosta-Humanez P. B. and Blazquez-Sanz D., Discr. and Contin. Dynam. Sys. Series B, 10, 2&3, 265-293 (2008). 

[3] Stachowiak T., Szydlowski M. and Maciejewski A. J., J. Math. Phys. 47 032502 (2006). 

[4] Cooper P., Khare A. and Sukhatme U., Phys. Repts. 251 267-385 (1995). 

[5] Kaplansky I., "An Introduction to differential algebra", Hermann, Paris (1957). 

[6] Kovacic J., J. Symbolic Comput. 2(1): 3-34 (1986). 




(13) 




(14) 



IV. ACKNOWLEDGEMENTS 



